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ABSTRACT 

We present results of a search for a young stellar moving group associated 
with the star HD 141569, a nearby, isolated Herbig AeBe primary member of a 
5±3 Myr-old triple star system on the outskirts of the Sco-Cen complex. Our 
spectroscopic survey identified a population of 21 Li-rich, <30 Myr-old stars 
within 30° of HD 141569 which possess similar proper motions with the star. 
The spatial distribution of these Li-rich stars, however, is not suggestive of a 
moving group associated with the HD 141569 triplet, but rather this sample 
appears cospatial with Upper Scorpius and Upper Centaurus Lupus. We apply 
a modified moving cluster parallax method to compare the kinematics of these 
youthful stars with Upper Scorpius and Upper Centaurus Lupus. Eight new 
potential members of Upper Scorpius and five new potential members of Upper 
Centaurus Lupus are identified. A substantial moving group with an identifiable 
nucleus within 15°(~30 pc) of HD 141569 is not found in this sample. Evidently, 
the HD 141569 system formed ~5 Myr ago in relative isolation, tens of parsecs 
away from the recent sites of star formation in the Ophiucus-Scorpius-Centaurus 
region. 
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INTRODUCTION 



The HD 141569 stellar system initially garnered the interest of iRossiterl (119431 ) as a 
potential t riple star system. The ex istence of this group was confirmed over half a century 
later when [Weinberger et al.l (120001 ) showed that the stars' position angles relative to one 
another had not changed since Rossiter's 1938 observations, thus confirming common proper 
motions. They also noted that all three stars are consistent with being the same age, indi- 
cating a comoving and coeval system. This triumvirate is particularly i nteresting because 



all co mponents are quite young; the system was found to be 5±3 Myr (jWeinberger et al. 



2000l ). and recently the A st ar's age was const rained using surface gravities and effective 
temperatures to be ~4.7 Myr (IMerm et al.ll2004j ). HD 141569 itself is near enough to resolve 
its large, dusty disk, and near-infrared signatures i ndicate perturbations t o the disk's struc- 
ture which could be explained by planet formation (jWeinberger et al.lll999l ). Were a coherent 
group present, it would therefore be distinctly possible to similarly observe additional young, 
disk-bearing stars. Furthermore, surveying a coeval sample at this distance would be useful 
for determining disk frequency and mechanisms by which they arise and evolve. 

In some cases, seemingly isolated young stars have later been found to have an entourage 
of other low-mass yo ung stars which const itut e a st ellar association (e.g. HD 104237, HR 
4796, (3 Pictoris; cf. iMamajek et al.l l2002l : iLiI |2005| ). The lower mass members of these 
associations have been useful, indeed critical, for estimating the age of the massive star 
in question, as often the massive star is on the main sequence and uncertainty in its HR 
diagram position is substantial enough to impede an isochronal age estimate. Any additional 
association members should be identifiable by their similar space motions. The youth of these 
objects makes them useful for studies of circumstellar disks and stellar evolution. In this 
region of sky rich with already identified associations such as p Ophiucus, Upper Scorpius 
(US) and Upper Centaurus Lupus (UCL), placing an "HD 141569 Association" into context 
with surrounding moving groups could aid in the understanding of star formation histories 
in giant molecular clouds. 

So motivated, in this paper we seek to identify a sar nple o f stars associated with 



HD 141569. Previous studies (e.g. iMamajek et al.l l2002l : iLiI |2005| ) have shown that cata- 



log searches based on X-ray activity, proper motions, and distance criteria can yield new 
members of associations. Here we also identify candidate members with a catalog search for 
X-ray active stars in the vicinity of HD 141569 possessing proper motions consistent with 
HD 141569 ( §2.ip . From new spectroscopic observations ( §2.2p . in fj3]we identify a subset 
of 21 stars which we claim as youthful using Li I equivalent widths; furthermore, we also 
note the presence of Ha emission as an interesting quantity usually indicative of chromo- 
spheric activity. Lacking direct distance measurements for these stars, we derive distances 
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using a modified vers i on of the tradit ional moving cluster parallax method (Ide Bruijndll999 
de Zeeuw et al.lll999l : lMamajekll2005l ). with which we further determine the probability that 
the sample stars are comoving with HD 141569 ( §3 ■41) . 

Placing these stars on a Hertzsprung-Russell diagram (hereafter HRD or HR diagram, 
§4.3p . we find that the stars' distances — and hence ages — are most consistent with their 
strong lithium abundances when we derive their moving cluster parallaxes assuming co- 
movement with Upper Scorpius or UCL associations of young stars. In [IS] we summarize our 
study and present thirteen newly identified Sco-Cen members. 



2. DATA 



2.1. Target Selection 



Potential targets for observation were found by two catalog queries for stellar properties 
indicative of youth and comovement with HD 141569. In the first sea rch, 10 stars were 
found through a query of the Hipparcos catalog (jPerryman fc ESAl 19971) for objects within 
10° of HD 141569 with similar distances (99±8pc, as reported in IPerrvman fc ESAlll997h 
and proper motions. Eight objects from the ROSAT Faint Source Catalog (jVoges et al. 



20001 ) supplemented this sample, as previous low-resolution observations (A. Weinberger, 
unpublished data) showed evidence for spectroscopic signatures indica tive of youth. A se cond 
search was performed in which the ROSAT Bright Source Catalog (jVoges et al.lll999l ) was 
probed for X-ray sources w ithin a 30° radiu s of HD 141569. Of the 1,114 resulting targets, 
~400 sources had Tycho-2 (lH0g et al.ll2000l ) catalogued proper motions. We required proper 
motions to be within ±15 and ±20 mas yr~^ in RA and Dec of HD 141569's proper motion, 
respectively, based on the range of pro per motions observed in the widely disp ersed and 
similarly aged TW Hydrae Association (IZuckerman et al.ll200ll : IWebb et al.lll999l ). Finally, 
we required that the ROSAT sources not be extended and that they be cospatial with 
the Tycho coordinates. This procedure resulted in ~70 objects desirable for further study 
to establish youth and space motions. A map of the observed sample is shown in Fig. [U 
members of other nearby associations are also shown for added context and perspective. 



2.2. Observations 

Forty-nine (49) stars of our input catalog were observed spectroscopically over the course 
of three observing runs in order to measure spectral types, radial velocities, and Li I and Ha 
equivalent widths. Four targets were found to be close visual binaries so their companions 
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were observed as well; of these four, two were determined to be background objects on the 
basis of apparent magnitudes and spectral type in comparison to the primary and are thus 
not discussed further. An observing log of the 51 stars (49 targets + 2 visual companions) 
which form our final sample for this work is presented in Table [TJ Photometry, proper 
motions, and parallax measurements from the literature are documented in Table [2l 

On 2001 June 18 (UT), the ten Hipparcos selected targets and five FSC targets were ob- 
served using the Hamilton Echelle Spectrometer at Lick Observatory. The Hamilton echelle 
covers a wavelength range of ~3500-10000A with a resultant resolving power of 60,000 at 
6OOOA when using a slit width of 1.2". During April 2002, a further 36 target stars of the 
second catalog search were observed using the echelle spectrograph on the Irenee du Pont 
telescope at Las Campanas Observatory. The du Pont spectrograph observes a wavelength 
range of ~3700-9800A with a resolving power of 40,000 at a slit width of 0.75". On both 
observing runs, comparison arc spectra were taken to establish a pixel to wavelength cali- 
bration. At LCO, Thorium-Argon lamp spectra were taken between each object exposure 
because the spectrograph is situated at the Cassegrain focus, whereas the Coude-fed Hamil- 
ton echelle required less frequent arc calibration, with arc spectra be taken at the beginning 
and end of the night. 



2.3. Reduction Procedure 

The echelle data were reduced using to perform standard spectral reduction 

procedures. Instrumental effects were accounted for via measurement of bias level and read 
noise, variations in pixel-to-pixel CCD response were removed in the flat fielding process, 
and dead columns were identified (which then prevented measurement of the Ha 6563 A 
line). A basic reduction process was then followed to locate and extract the echelle orders 
as well as remove scattered light. The Thorium- Argon arcs were similarly extracted, their 
features identified so as to wavelength calibrate the object spectra. Examples of IRAF 
routines employed include apall, ecidentify, refspec, and dispcor. 



""^Image Reduction and Analysis Facility, IRAF, is distributed by the National Optical Astronomy Ob- 
servatories, which are operated by the Association of Universities for Research in Astronomy, Inc., under 
cooperative agreement with the National Science Foundation. 
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3. ANALYSIS 

We aim to identify the subset of our sample which shows evidence for youth consistent 
with membership to a young, ~5 Myr-old association. Additionally, we wish to assess 
kinematic properties and test for comovement. To these ends, from the literature, we obtain 
effective temperatures as well as near infrared photometry and proper motions. From the 
spectra, we measure Li I equivalent widths as well as radial velocities. Utilizing these data, 
we estimate ages and select a youthful, high lithium sample, and we test for kinematic 
similarity against the velocity models of nearby, young moving groups. 



3.1. Effective Temperatures 



Spectral types for the ma.iority of our samp l e are reported in t he hterature (lHouklll982 
Houk fc Smith-Moorelll988l : iHouk fc Swiftlll999l : iTorres et al.ll2006l). These are converted t o 
Teff using the main-sequence spectral type-Teg relationship of iKenyon fc Hartmanru (119951 ). 
We supplement these with spectral types determined from low-resolution spectra obtained 
with the Lick KAST spectrograph (A. Weinberger, unpubli s hed data) as well a s T^s deter- 
mined from the Fe I to Sc I line ratio (cf. IStassun et al.ll2004l : ISteffen et al.ll200ll ) observed in 
our high- resolution spectra ( §2.2p . The line ra tios measured for our sp ectral type standards 
are in good agreement with the calibration of iBasri fc Batalhal (Il990l ): we thus adopt their 
line ratio-spectral type scale in assigning types to our sample. In some cases, primarily for 
more massive stars, the Fe I and Sc I lines were not present or could not be measured with 
confidence above the noise level. T^s values are summarized in Table [31 Where multiple T^s 
values are available, they generally show good consistency with one another to within ~ 300 
K (corresponding to ~ 2 spectral subclasses). For our final set of effective temperatures, 
we adopt literature spectral types where available, line ratio spectral types if not. For one 
object, HD 1573 10B, neither were ava ilable, and thus we interpolated its 2MASS (2 Micron 



All Sky Surve y. ISkrutskie et al.l 120061) (H — Ks) color over the effective temperature-color 



relationship of Kenyon fc HartmannI jl995 ). The set of adopted temperatures is reported in 



the final column of Table [31 

In Fig. [21 we show these effective temperatures as a function of th e objects' observed 
2MASS (H—Ks) colors; also plotted is the Tef[ — {H—Ks) relationship from lKenyon fc Hartmann 



(Il995l ). For comparison, the standard stars observed (Table [T|) are also shown in this pa- 
rameter space. The observed {H — Kg) colors follow the expected relationship with T^s with 
a scatter of ~ 300 K, consistent with the scatter in Teg from spectral types above. This 
indicates that our sample in ge neral suffers relativel y little extinction. Indeed, radio survey 
column-density measurements (iKalberla et al.ll2005l ) in the direction of our targets indicate 
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that E{H — Kg) reddening toward our sample should be < 0.1 mag. Via comparison with 
expected intrinsic {H — Kg) for a star of given effective temperature, we derive Kg band 
extinction values, A^s^ and deredden the sample accordingly. {H — Kg) colors are plotted in 
Fig. [31 lines illustrate dereddening by connecting colo rs before and after der eddening. Also 
displayed for comparison are the dwarf sequence from lBessell fc BrettI (119881 ) as well as red- 
dening vectors assuming the standard ratio of total-to-selective extinction Ry = 3.12. For 
visual clarity, we only display the 21 targets identifed as lithium rich ( §3.2p . The X^-band 
extinction corresponding to the applied dereddenings is in all cases A^s < 0.15 mag, with 
the exceptions of 2MASS J17215666-2010498 and TYC 6191-0552-1 (objects #19 and #2 
in the data tables) for which A^s = 0.33 mag and Aks = 0.16 mag, respectively. We have 
checked Spitzer 24/im data (A. Weinberger, private communication) and find that while 
most of the high Lithium sample lacks 24/im excess, 2MASS J17215666- 2010498 has a sub- 
stantial excess. TYC 6191-0552-1 was analyzed by iMeyer et al.l (120081 ) and was found to 
have moderate 24/im excess. For both objects, apparent excess in the H and Ks bands in 
concert with 24um excess confirms disk presence and thus we cannot and do not apply a 
standard interstellar dereddening law. 



3.2. Lithium Equivalent Width 



We measured the equivalent width (EW) of the A6707 line of Li I from our spectra 
using the IRAF routine Splot. For each star, EW measurements were obtained by both 
directly integrating the fiux in the line and by calculating the area of a best-fitting Gaussian. 
Our measured EW values include contributions from the small Fe I+CN line at 6707. 44A, 
leading to measured Li I EWs that are rep resentative of a sl i ghtly (10-20 mA) over-estimated 
photospheric Li presence. For instance ISoderblom et al.l (119931 ) report that this Fe line 
blend has an EW = [20(B-V)o-3]mA for main sequence solar-type stars. We correct for 
contamination following this prescription and find the median value for the sample is 10 
mA; we report in Table [3] the Li I EWs for each target via both measurement methods as 
well as the Fe line blend contribution. For all targets, the RMS of the difference between the 
EWs determined using both methods is 18.4 mA. With the aim of conservatively selecting 
a sample, we utilize the lower of the two measured values; these are plotted in Fig. |H 

To identify the young stars in the sample, Li I EWs were compared to the upper envelope 
of EWs as a f unction of T^ff reported in t he literature for the ~ 30 Myr-old clusters IC 2602 
and IC 2391 (IRandich et al.l 119971 . 1200 ll . our Fig. H]). Twelve stars are found to have Li I 
EWs above this envelope, compelling evidence of youth. In what follows we refer to these 
12 stars as the "High Li" sample. 
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Also of interest, the Ha profiles of several objects are in emission in our spectra and 
some possess double-peaked profiles. Another 6 stars show elevated lithium levels (> 200mA 
but are below the threshold shown in Fig. H]), placing them in the upper envelope of the IC 
2602 and IC 2391 loci. An additional two stars have temperatures greater than 7500K and 
equivalent widths above the locus; although these are potentially older stars which simply 
lack deep enough convective zones to deplete primordial lithium abundances, we include 
them in the analyses for completeness. We also include star ^13 in this group as it is in a 
double system separated by ~1.2"while the seeing that night was ~1.5", making Li line filling 
likely. In Fig. HI this star is plotted with its measured equivalent width doubled (denoted 
"13b") to demonstrate which sample group it would potentially belong to. In what follows, 
we refer to these 9 stars as the "Moderate Li" sample, two of which show double-peaked Ha 
in emission (Table [3]). 

These 21 stars, which are the most likely in our sample to be of comparable age to HD 
141569, will be the focus of the remainder of our analyses. For ease in tracking these stars 
through our analysis, they are labeled with a running numerical identifier in the data tables 
and figures, and all EW and line profile information is reported in Table [3l 

3.3. Radial Velocity 

Heliocentric radial velocities were obtained using the IRAF task fxcor to cross correlate 
each target spectrum against the radial- velocity standard star of closest spectral type (Table 
[1]). Four echelle orders spanning the wavelength ranges 6025-6150A, 6150-6275A, 6625- 
6750A, and 5120-5220A were employed as they contain many deep metallic lines and little 
or no telluric contamination. In Table H] we report the mean radial velocities from the four 
orders. 

To determine the extent to which our radial-velocity measurements may be affected 
by the vsini and signal-to-noise (S/N) of our target spectra, we performed a Monte Carlo 
simulation in which a narrow-lined, high S/N standard star spectrum was randomly degraded 
one hundred times. This process created artificially noisy spectra at S/N levels of 10-100, 
well representing the full range of S/N found in our sample. Each degraded spectrum was 
artificially broadened and cross-correlated against its original high S/N spectrum. These 
degraded spectra were furthermore cross-correlated against the other standard stars to assess 
the effects of spectral-type mismatch on the resulting radial velocities. We find that these 
effects are negligible (i.e. affecting the resulting radial velocities by < 1 km s~^) unless 
vsini > 70 km s~^ or S/N < 30. As all of our target spectra have S/N > 30, only very fast 
rotators are potentially affected (by up to 2.3 km s~^ for vsini = 100 km s~^). Measured 
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V sin i and radial velocity values are documented in Table HI In addition to these effects, we 
note that on an aperture-to-aperture basis, errors in wavelength calibration could affect the 
measured radial velocity; this can be quantified as the standard deviation of the mean radial 
velocity measured from the four selected apertures. The final radial-velocity uncertainties 
quoted in Table H] are the quadrature sum of the internal uncertainty (the standard deviation 
of the four spectral orders used) and the uncertainty arising from rotational broadening. 



3.4. Moving Cluster Parallaxes 



With observed proper motions and measured radial velocities for the 21 "High" and 
"Moderate" lithium stars, a kinematic picture of the sample is almost complete. Tycho-2 
proper motions were used for consistency throughout (Table |2]), the only exception be- 
i ng one of the two clo se visual binary companions, HD 157310B, for which only UCAC2 
( IZacharias et al.l |2004J ) proper motions were available. Proper motion data are not avail- 
able for the other close companion, 2MASS J17215666-2010498; in what follows, we assume 
common proper motion with its primary star, TYC 6242-0104-1. Hipparcos pa rallaxes are 
unavailable for the high lithium stars, thus a moving cluster parallax method (Ide Bruijne 
19991 ) provides a means for determining their parallaxes and hence their distances. With 
distances, it can be tested then whether these objects are consistent with being members of 
a coherent moving group. 



Our procedure is rooted in the derivation of Ide Bruijnd (119991 . see their §2 and refer- 
ences therein). The process is executed assuming a velocity vector, and hence a convergent 
point, for the moving group to which we are testing membership. In analyzing the spatial 
distribution of the 21 stars of our youthful sample, we note they are all in closer proximity to 
Upper Sco and UCL than HD 141569 (see Fig. [5]). The youthful sample is highly spatially 
separated from HD 141569, and these separations indicate two kinematic issues. First, it 
is unlikely for objects with such large separations to be comoving. Second, had these stars 
indeed formed together, large initial velocities (~6 km s~^, inconsistent with the observed 
1-2 km s~^ velocity dispersions of young associations) would be required to bring about the 
separations presently observed after ~5 Myr of motion. As it is unlikely these objects are 
associated with HD 141569, we require estimates of the mean velocity vectors for Upper 
Sco and UCL. The velocity vector for Upper Sco is adopted from lMamajekl (120081 ): UVW = 
[—5.2, — 16.6,— 7.3] km s~^. This vector incorporates a mean radial velocity for 120 Upper Sco 
members, an improvement over prior velocity vectors which solely rely upon proper motion 
and parallax information. For UCL, the velocity model used for comparison is derived from 



the median position, proper motion, and radial velocities of UCL members ( Ide Zeeuw et al. 
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1999h : UVW = [-5.4,-19.7,-4.4] km s'^. We calculate these UVW vectors for Upper 
Sco and UCL to precisions of ~ ±0.3 and ~ ±0.4 km/s respectively, but note, there ex- 
ist discrepancies between UVW vectors derived by various authorso The reason for the 
systematic differences between published convergent points and velocity vectors for the OB 
subgroups is not completely clear. The leading candidates for these systematic differences 
are unaccounted-for expansion of the subgroups, and the probable presence of unresolved 
spatial and kinematic substructure within the subgroups. For robustness, we include all 
available radial velocity measurements in derivation of UVW vectors. 

For each of the 21 stars in our "High Lithium" and "Moderate Lithium" samples, we 
derive moving cluster parallaxes using the HD 141569 and Upper Sco velocity vectors (Table 
HJ. The formalism for this is: 

^= ■ /XX (1) 

V sm(Aj 

where zu is the parallax, A is 4.74 km yr s~^ (the ratio of one AU in km to a Julian Year 
in s), /Lt„ is the parallel component of proper motion (proper motion in the direction of the 
convergent point), v is the velocity of the group in km s~^, arid A is the angular separation 



between the star and the convergent point (formula 1 . lMamajekll2005l ) . An additionally useful 
parameter, the comovement probability, can also be calculated. Comovement probability is 
defined as 1 — Pi, where P± is the likelihood that the star's proper motion is entirely 
perpendicular to the direction of the convergent point; the projection of proper motion in 
this direction is denoted Hr, and a fi^ close to is indicative of comovement. 



4. RESULTS 

4.1. Distances, comovement probabilities, and membership 

The spatial proximity of our youthful sample stars to Upper Sco and UCL (Fig. 
suggests that these objects are not likely related kinematically to the farther away HD 141569 
system. It is important to stress at this juncture that derived comovement probabilities are 
not absolute probabilities per se; their derivation depends directly on the velocity model 
assumed a priori. We do know with certainty that the stars in our sample are young (by 
virtue of their high Li abundances) and that they are moreover in projected proximity to 
other stars known to be young, nearby, and comoving (Fig. [1]). Thus there is a strong "prior" 



^For example , 



Madsen et al 



(120021 ) calculate for Upper Sco UVW = [-0.9,-16.9,-5.3] km s-\ 
de Bruiine et al.l (j200l[ ) cite UV W = [4.1,— 17.9 , — 3.7] km (with model-observation discrepancy param- 



eter "g" set to equal nine) while Ide Zeeuw et al.l (|l999l ) report UVW = [0.0,-16.1,-4.6] km s" 
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favoring the velocity models that we have chosen to test. Still, the comovement probabilities 
reported in Table H] should be regarded as measures of consistency with the assumed velocity 
models, not proof of membership. We therefore adopt the very simple criterion of spatial 
proximity to a group in application of velocity modeling, and report the resulting distances 
and comovement probabilities for objects when tested against the velocity vectors of Upper 
Sco and UCL. Parallax distances and comovement probabilities calculated as previously 
described ( §3.41) are reported in Table HI 

Based on two simple criteria, youth determined via measurement of the A6707 line and 
spatial positio n, objects 2, 6, 7 , 8, 9, 11, 13, 15, and 21 lie within the Upper Scorpius "box" 



as defined by Ide Zeeuw et al.l (119991 ). Similarly, stars 4, 5, 14, 17, and 20 appear to be 
UCL members. These "spatial matches" are summarized in column 3 of Table O Due to 
the similarities of velocity vectors in the Ophiucus-Sco-Cen region, it is unsurprising that 
in some cases, objects we deem US or UCL members have higher comovement probabilities 
when tested against the velocity vector of the other group. Factors which create blurring of 
kinematic boundaries include internal velocity dispersions inherent to a given moving group 
and observational uncertainties which then propagate into the convergent point solution. We 
thus take {£, b) position as the strongest indicator of group membership and then examine 
comovement probabilities as a supplement. Outside of the Upper Sco box, stars 3, 12, and 
18 have high comovement probabilities with Upper Sco. We would present 12 and 18 with 
some caution as Upper Sco members, as they are within a few degrees of the most extended, 
already known Upper Sco members. 

Object 3 is almost ten degrees away from the southernmost US stars and thus its 
association with Upper Sco is also dubious. For these three objects, we tentatively suggest 
Upper Sco membership and denote their membership in Table [5] as "US?" In two cases 
we note objects with low comovement probabilities with their spatially matched groups: 
objects 1 and 19 do not have velocities consistent with US and are spatially inconsistent 
with being UCL members. Object 15, while spatially coincident with Upper Sco, has low 
enough comovement probability to be suspect. Finally, stars 10 and 16 have high comovement 
probabilities with Upper Sco, but appear to be too far away in (£, b) space to be considered 
part of Upper Sco. These remaining five objects we also classify as being of "Indeterminate" 
membership. In summary, the total number of new Upper Sco members presented here is 
eight, and five new members of UCL are also identified. 
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4.2. Space Motions 

To illu strate kinematic associa t ion in a familiar way, we could utilize the transformation 



matrices of iJohnson fc SoderblomI (119871 ) to calculate UVW space motions for the sample 
stars. UVW motions, however, depend on distance, a quantity we have obtained via assump- 
tion of comovement with a given UVW vector. The resulting UVW plot is thus degenerate 
and does not provide additional criteria by which we can further examine association. 

As an additional check on the application of each velocity model, given an assumed 
velocity vector, radial velocities for each object can be predicted based on their proper 
motions and positions. We find good consistency between the predicted radial velocities 
and those measured when comparing measured radial velocity to predicted radial velocity 
for whichever velocity model we would naively expect need be applied given simple spatial 
proximity to a given moving group. Illustrating the radial velocity structure of the sample in 
context with nearby groups can be a measurement-based, assumption-free way of analyzing 
space motions. In the selection criteria we constrained proper motions to agree with those of 
HD 141569 within a wide range of values which includes proper motions generally observed 
in Upper Sco. Measured radial velocities and projected radial velocities for the Upper Sco 
velocity vector are plotted as a function of galactic longitude in Fig. [61 Most notably, the 
entire "High Li" sample agrees well with the predicted radial velocities of Upper Sco within 
~2-3cr. The farthest outlying points are from the "Moderate Li" sample. 



4.3. H-R Diagram 



In Fig. [7] we show the placement of the sample stars on three HR diagrams to illustrate 
shifts in M^^ magnitude due to changes in distance. Absolute magnitudes were calculated 
from the observed 2MASS K5 magnitudes (Table [2]) and one of three distances. Uncertainties 
in M-Ks are the propagated errors of the 2MASS photometry together with the formal errors 
in distance. The uncertainty in T^s is taken to be two spectral subtypes (see §3.11) . To derive 
ages for our sample stars we also show the pre-n i ain-se quence (PMS) evolutionary tracks of 
Baraffe et all Jl998h and IP'Antona fc MazziteUil Jl997l ). 



In the upper panel, we illustrate placement of the sample on the HRD when we apply the 
distance to the HD 141569 system to every individual object. The isochronal ages inferred 
for most stars in our sample using the HD 141569 mean distance are in general older (~ 30- 
100 Myr) than what would be expected for the stars based on their lithium abundances and 
comovement with the HD 141569 system (age 5±3 Myr). In concert with the lack of spatial 
proximity, we further rule out the potential for a coeval, coherent moving group near HD 
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141569. 



In contrast, the inferred ages using the Upper Sco mean distance (~145pc, effectively 
equivalent that of UCL, ~142pc) are entirely consistent with the expectation of <30 Myr as 
imposed by the Li EW measurements. All stars appear on or above the 30 Myr isochrone, 
save objects t^IO and #16, which, despite their high comovement probabilities with Upper 
Sco, do not appear to be in close enough spatial proximity to be members of that moving 
group. In the third HR diagram, applied distances are determined by which velocity model, 
US or UCL, provides the highest comovement probability with a given object. This particular 
representation is not only mostly consistent with the age range expected from Li I presence, 
it also "correctly" places the higher mass objects closer to the ZAMS; particularly, objects 
10, 16, and 17 have derived distances which make them appear to be ZAMS stars rather 
than anomalous objects far above or below the theoretical isochrones. 

Scatter in an HRD generally can be attributed to many factors; the radial extent of a 
young association (jMamajekll2005l . e.g., TW Hydra ~55pc), observational errors, or even the 
choice of evolutionary tracks can gener ate shifts and enhance spread in isochronal age of a 
sample expected to be coeval. Using the lD'Antona fc Mazzitellil (119971 ) evolutionary models, 
our sample stars all appear to have isochronal ages <10 Myr. In general, these track s appear 
shifted by ~400K to higher temperatures with respect to the iBaraffe et al.l (119981) tracks. 
For an illustration of these track-based discrepancies, see Simon et al. (j2000 ). In spite of 
obstacles posed by apparently discrepant HR diagrams, we can say with confidence based on 
Li I presence that these objects are indeed young, <30 Myr. The HR diagram, when applying 
high comovement probability derived distances, provides a higher degree of confidence in 
adopting these distances as the ages are indeed as expected from Li measurements. 



4.4. Is HD 141569 Related to Upper Sco? 

HD 141569 is an apparently isolated system located within tens of degrees (and parsecs) 
of known sites of recent (<5 Myr) and ongoing star formation, all apparently associated with 
the Sco-Cen star forming complex {d = 100-200 pc; Prei bisch fc Mamajek, 2008, in press). 



which appears contiguous with the Aquila Rift regions (iDame et al.l 119871 . i ~ 30°;). HD 
141569 appears to agree with the projected velocity model of Upper Sco (see Fig. [6l but 
we can rule out the possibility of HD 141569 originating or being kinematically associated 
with Upper Sco. The hypothesis that HD 141569 could have been ejected at high velocity 
from a known high density stellar nursery can be strongly discounted on two grounds. First, 
HD 141569 has two low-mass companions at wide separation (~10^ AU, with likely orbital 
motion of ~1 km s~^). A velocity kick of >2-3 km s~^ to either A or B-l-C would have 
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likely disintegrated the system. A low velocity ejection (<2-3 km s~^) would have placed the 
birth site within <10-15 pc (<5°-7°), but no such known young clusters or molecular clouds 
appear there. 

Position and velocity information for these nearby groups was entered into an orbit 
code which employs the epicyclic approximation. The separations between these groups and 
HD 141569 were evaluated during the past 10 Myr; combined distance and velocity vector 
uncertainties result in <15 pc uncertainties over this time frame. Presently, HD 141569 is 
~55 pc away from the center of Upper Sco, and was only slightly closer at its minimum 
separation of ~53 pc (~2.7 Myr ago). In UVW, the only substantial difference is in the 
W component of velocity: w hile Upper Sco has negligible vertical motion with respect to 



the Local Standard of Rest (iMamajekl |2008| ). HD 141569 is moving northward out of the 



disk at ~5 km s~^. This anomalous W component of motion is discrepant with any known 
molecular cloud or star forming region near HD 141569. Further, given the kinematic data 
and the isochronal age of HD 141569 (~4-5 Myr), it appears that HD 141569 could not 
have formed from any of the known sites of recent star-formation in its vicinity. The list 
of excluded birth-sites includes Upper Sco, UCL, Lower Centaurus Crux (LCC), and the 
Ophiuchus, Corona Australis, and Lupus clouds. 

Combining the kinematic, position, and age data, we conclude that the HD 141569 
triple system likely formed in isolation or with a small entourage of companions in a cloud. 
It appears to have formed ~25 pc closer to the Galactic plane than its present position, 
and its anomalously large W velocity component has carried it to its modern high latitude 
position (6 ~ +37°, ~90 pc above the Galactic plane). 



5. Summary and Conclusions 

We have identified a group of 21 PMS stars within 30° of HD 141569 on the basis of 
strong Li I absorption and, in 9 of those 21 cases. Ha in emission. These stars were selected 
through a joint catalog search for X-ray sources with spatial and proper-motion characteris- 
tics similar to those of HD 141569, a B9.5Ve star at 116 pc that harbors a circumstellar disk 



and for which two low-mass companions had previously been identified ( IWeinberger et al. 



20001 ). For these 21 stars, we have applied a moving cluster parallax technique to proper- 



motion data from the literature. 

Table [5] outlines our final membership assessments: we present eight potential new 
members of Upper Sco and five new potential members of UCL. These stars possess lithium 
presence consistent with youth and futhermore appear youthful on the HR diagram. Pri- 
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marily we utilize spatial position as the principal criterion for determining membership and 
supplement that with comovement probabilties from the moving cluster parallax derivation. 
Additionally, we examine the motion of the HD 141569 system away from the galactic mid- 
plane over its lifetime. The system, surprisingly, appears to have formed in isolation, well 
outside of presently known star forming regions and molecular clouds. 
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Table 1. Observing Log 



Object Name 


Right Ascension 


Declination 


Observation Time 


Integration 


s/n(i) 


Comment(s) 




[J2000] 


[J2000] 


[UT] 


Time [s] 






2001-06-18 : UCO / Lick 


Alpha Boo 


14:15:39.67 


+19:10:56.7 


04:17:13.0 


1 


200 


K2III Vr Standard 


HD 137396 


15:26:05.91 


-11:41:55.7 


04:28:35.0 


720 


127 




RHS 48 


15:23:46.0 


-00:44:25 


04:52:38.0 


1500 


158 




HD 138969 


15:35:47.41 


-12:51:32.9 


05:25:23.0 


720 


105 




xlU 14Uo / 4 


io;44;zo.oU 


-Oo:oO:lo.o 


r\r^, A coo n 


Ton 


132 




HD 141612 


15:50:11.35 


-06:32:14.6 


06:05:08.0 


1500 


178 




HD 141693 


15:50:11.86 


-1-05:57:17.2 


06:37:21.0 


360 


105 




HD 142987 


15:58:20.558 


-18:37:25.09 


06:50:42.0 


1800 


75 




TYC 6191-0552 


15:58:47.724 


-17:57:59.79 


07:30:24.0 


2400 


31 




HD 143332 


16:00:02.65 


-07:07:21.8 


08:17:18.0 


600 


126 




HIP 79354 


16:11:44.57 


-1-00:47:38.6 


08:34:18.0 


1000 


97 




HD 145551 


16:11:56.93 


-10:37:20.3 


08:58:01.0 


600 


98 




BD -06 4414 


16:22:17.866 


-06:23:03.93 


09:15:34.0 


600 


72 




HD 143810 


16:02:47.47 


-10:37:56.7 


09:32:25.0 


600 


72 




HD 144726 


16:07:44.91 


-11:55:13.1 


09:49:16.0 


600 


145 




HD 145169 


16:09:50.80 


-11:07:56.5 


10:06:16.0 


450 


95 




HD 177178 


19:03:32.25 


+01:49:07.6 


10:20:40.0 


300 


175 


A4V Vr Standard 


HD 187691 


19:51:01.64 


+10:24:56.6 


10:33:26.0 


240 


184 


F8V Vr Standard 


16 Cyg/HD 186427 


19:40:32.06 


+50:31:03.1 


10:45:50.0 


240 


197 


G2V Vr Standard 


HD 184467 


19:31:07.97 


+58:35:09.6 


10:59:00.0 


300 


139 


K2V Vr Standard 


2002-04-17 ; LCO 


HD 80170 


09:16:57.08 


-39:24:00.5 


01:07:53.4 


30 


107 


K5III Vr Standard 


HD 102870 


11:50:41.72 


+01:45:53.0 


01:16:31.1 


20 


131 


F9V Vr Standard 


TYC 6141-0525-1 


14:06:06.475 


-18:10:37.41 


03:17:19.9 


1080 


81 




TYC 0909-0125-1 


14:19:33.692 


+08:38:09.59 


03:43:18.4 


1080 


79 




TYC 7312-0236-1 


15:12:44.473 


-31:16:48.10 


00:00:00.0 


1500 


89 




TYC 7327-0689-1 


15:37:51.344 


-30:45:16.15 


04:37:09.7 


1500 


77 




TYC 5003-0138-1 


15:26:52.727 


-00:53:11.74 


05:09:45.3 


1500 


112 




TYC 0937-0754-1 


15:38:00.94 


+11:34:57.30 


05:41:42.1 


1500 


75 




TYC 6781-0415-1 


15:41:31.212 


-25:20:36.39 


06:16:54.2 


900 


86 




TYC 6790-1227-1 


15:48:02.922 


-29:08:37.02 


06:40:12.4 


1500 


101 




HD 142016 


15:-53:12.969 


-30:46:44.17 


07:26:45.0 


255 


58 




TYC 0376-0769-1 


16:19:59.977 


+04:36:46.33 


07:42:44.8 


1200 


85 





Table 1 — Continued 



Object Name 


Right Ascension 


Declination 


Observation Time 


Integration 


S/N(i) 


Comment(s) 




[J2000] 


[J2000] 


[UT] 


Time [s] 






TYC 7346-1182-1 


16:42:07.711 


-30:38:37.87 


08:12:47.3 


1200 


117 




HD 153439 


17:00:42.963 


-27:25:17.55 


08:37:23.9 


750 


136 




CD -25 11942 


17:06:01.190 


-25:20:30.40 


08:57:16.9 


900 


150 




TYC 6817-1757-1 


16:49:35.996 


-27:28:07.78 


09:18:32.5 


1350 


80 




HD 188376 


19:55:50.36 


-26:17:58.2 


09:47:46.4 


65 


109 


G5V Vr Standard 


HD 165341 


18:05:27.29 


+02:30:00.4 


09:58:15.0 


45 


206 


KOV Spectral Type Standard 


HD 209290 


22:02:10.27 


01:21:00.8 


10:06:18.1 


750 


52 


]\10.5\'^ S])cctral Typo Standard 


2002-1)1-18 : LCX) 


TYC 5022-0263-1 


15:48:40.942 


-03:10:14.22 


06:22:52.4 


1500 


116 




HD 143358 


16:01:07.927 


-32:54:52.53 


06:54:12.0 


900 


117 




HD 144713 


16:08:05.223 


-24:55:33.2 


07:14:23.6 


900 


153 




HD 144732 


16:08:26.317 


-28:25:49.3 


07:34:59.7 


1200 


105 




TYC 6806-0888-1 


16:23:46.978 


-28:50:02.43 


07:59:41.7 


1800 


100 




TYC 6803-0897-1 


16:29:49.911 


-27:28:49.96 


08:35:31.7 


1800 


99 




HD 148982 


16:32:50.876 


-28:20:39.89 


09:11:17.0 


1200 


115 




HD 157310 


17:22:23.477 


+04:45:24.83 


09:38:10.9 


1200 


99 


Double star - brighter component 


BD +04 3405B 


17:22:23.5 


+04:45:13.3 


09:59:37.2 


1200 


60 


Double star - fainter component 


2002-04-19 : LCD 


HD 141813 


15:51:54.393 


-26:22:05.52 


05:41:33.5 


900 


120 




HD 148396 


16:28:51.327 


-30:16:55.26 


06:06:32.7 


1500 


119 


Catalogued Double star 


TYC 0976-1617-1 


17:01:57.543 


-1-07:33:32.45 


06:40:26.6 


2400 


100 


SB (2) 


HD 154922 


17:09:07.521 


-18:15:22.78 


07:27:02.2 


450 


107 




2MASS J17215666-2010498 


17:21:56.67 


-20:10:49.8 


08:03:30.1 


3600 


52 


Double star - brighter component 


TYC 6242-0104-1 


17:21:56.048 


-20:10:51.97 


09:09:59.6 


3600 


71 


Double star - fainter component 


2002-04-20 : LCO 


HD 109524 


12:35:33.55 


-34:52:54.9 


04:11:22.8 


750 


238 


K2V Vr Standard 


HIP 75685 


15:27:42.638 


-02:45:18.56 


05:43:37.4 


900 


122 




TYC 6234-1287-1 


17:26:56.545 


-16:31:34.85 


06:06:12.0 


1800 


118 




TYC 7334-0429-1 


16:04:30.557 


-32:07:28.7 


06:43:18.0 


1500 


110 




TYC 5668-0365-1 


17:40:55.47 


-12:16:27.75 


07:14:42.1 


2400 


134 




HD 144393 


16:06:06.470 


-12:18:15.22 


08:01:37.6 


360 


109 




TYC 6214^2384^1 


16:19:33.958 


-22:28:29.41 


08:59:30.1 


1800 


99 




TYC 6215-0184-1 


16:29:48.700 


-21:52:11.89 


09:35:25.5 


2400 


129 


SB (2) 


HD 177178 


19:03:32.25 


-1-01:49:07.6 


10:24:04.7 


150 


109 


A4V Vr Standard 





Table 1- 


— Continued 




Object Name Right Ascension 
[J2000] 


Declination 
[J2000] 


Observation Time 
[UT] 


Integration S/N'^' Comment(s) 
Time [s] 





Note. — 

^ Approximated using Splot at ~6500Aand ~6700A. 

^■Suspected spectroscopic binary based upon broadened troughs of spectral features. 
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Tablc 2. Stellar Parameters 



Plot 


Object Name 




ij(i) 




Ma 




Parallax(3) 


ID 










[mas yr~^] 


[mas yr~^] 


[mas] 


A 


HD 141569 


6.872±0.027 


6.861±0.040 


6.281±0.026 


-18.3±1.1 


-20.5±1.1 


8.63±0.59 


1 


TYC 6242-0104-1 


9.963±0.027 


9.305±0.026 


9.151±0.024 


-11.7±3.6 


-13.7±4.0 




2 


TYC 6191-0552 


9.261±0.022 


8.535±0.042 


8.325±0.024 


-15.0±3.3 


-20.2±3.7 




3 


TYC 6234-1287-1 


8.659±0.025 


8.012±0.024 


7.829±0.020 


-10.1±2.6 


-39.0±2.7 




4 


TYC 7312-0236-1 


9.601±0.024 


9.079±0.024 


8.919±0.019 


-20.5±3.3 


-20.5±3.2 




5 


TYC 7327-0689-1 


9.300±0.024 


8.755±0.036 


8.563±0.019 


-19.7±3.2 


-27.5±3.1 




6 


TYC 6781-0415-1 


7.974±0.030 


7.367±0.033 


7.241±0.024 


-19.5±2.7 


-30.1±2.4 




7 


TYC 6803-0897-1 


9.275±0.024 


8.743±0.049 


8.648±0.025 


-15.6±2.5 


-28.6±2.5 




8 


TYC 6214-2384-1 


9.230±0.019 


8.659±0.036 


8.509±0.019 


-18.7±3.5 


-26.2±3.8 




9 


TYC 6806-0888-1 


9.216±0.025 


8.785±0.027 


8.659±0.026 


-13.4±3.0 


-27.5±2.7 




10 


BD -1-04 3405B 


9.759±0.022 


9.413±0.031 


9.262±0.019 


-5.3±1.5* 


-14.9itl.8* 




11 


HD 144713 


7.847±0.021 


7.538±0.034 


7.431±0.020 


-11.7±1.7 


-20.7±1.6 




12 


HD 153439 


8.073±0.020 


7.852±0.049 


7.729±0.047 


-6.7±1.6 


-28.6±1.6 




13 


HD 148396 


8.420±0.023 


8.095±0.019 


8.100±0.020 


-6.9±2.4 


-15.4±2.4 




14 


TYC 7334-0429-1 


9.168±0.018 


8.690±0.049 


8.565±0.021 


-17.5±2.2 


-25.5±2.2 




15 


TYC 6817-1757-1 


8.815±0.021 


8.350±0.042 


8.179±0.031 


-10.1±2.8 


-7.4±2.5 




16 


HD 157310 


9.160±0.022 


9.088±0.047 


9.006±0.021 


-5.4±1.5 


-12. Oil. 5 




17 


HD 142016 


6.785±0.020 


6.707±0.034 


6.622±0.018 


-26.4±1.2 


-38.8±1.3 




18 


CD -25 11942 


8.099±0.020 


7.661±0.029 


7.525±0.038 


-9.8±2.0 


-28.2±1.8 




19 


2MASS J17215666-2010498 


8.150±0.023 


7.187±0.047 


6.840±0.023 








20 


TYC 6790-1227-1 


9.212±0.023 


8.719±0.026 


8.624±0.023 


-20.4±2.8 


-26.0±2.3 




21 


TYC 7346-1182-1 


9.018±0.027 


8.663±0.053 


8.530±0.019 


-14.3±2.3 


-27.0±2.2 






HD 142987 


8.279±0.035 


7.774±0.063 


7.614±0.021 


-15.4it2.1 


-22.3±2.2 






HD 143358 


8.470±0.023 


8.164±0.036 


8.074±0.020 


-18.3±1.4 


-29.6±1.5 






BD -06 4414 


8.359±0.027 


8.043±0.038 


7.936±0.033 


-27.0±2.1 


-33.7±2.2 






HD 137396 


7.632±0.023 


7.493±0.033 


7.419±0.027 


-14.9±1.4 


-12.0±1.0 






HD 138969 


7.946±0.027 


7.693±0.040 


7.666±0.017 


-11.3±1.4 


-8.5±1.0 


10.69±1.4 




HD 140574 


7.623±0.018 


7.495±0.036 


7.459±0.029 


-27.2±1.3 


-30.2±0.9 


10.8±1.15 




HD 141612 


8.808±0.024 


8.497±0.057 


8.398±0.031 


-17.7±2.4 


-10.7±1.8 


9.62±1.74 




HD 141693 


6.886±0.020 


6.888±0.034 


6.828±0.023 


-29.7±1.2 


-26.5±0.7 


8.96±0.96 




HD 141813 


8.232±0.023 


7.963±0.036 


7.862±0.020 


-22.7±1.7 


-38.1±1.9 






HD 143332 


6.987±0.024 


6.783±0.044 


6.680±0.023 


-10.6±1.3 


-24.9±1.3 


9.96±1.36 




HD 143810 


8.953±0.030 


8.766±0.061 


8.643±0.019 


-24.8±1.5 


-13.7±1.5 


9.26±1.51 




HD 144393 


7.323±0.027 


7.099±0.046 


6.980±0.021 


-6.3±1.2 


-23.9±1.1 


10.97±1.14 




HD 144726 


7.527±0.027 


7.322±0.036 


7.244±0.026 


-16.7±1.4 


-19.4±1.3 


8.49±1.20 




HD 144732 


8.531±0.023 


8.231±0.051 


8.147±0.026 


-15.2±2.1 


-33.2±1.9 






HD 145169 


7.159±0.039 


6.912±0.040 


6.876±0.023 


-20.2itl.l 


-30.5±1.1 


12.75±1.78 




HD 145551 


7.906±0.029 


7.677±0.042 


7.600±0.018 


-20.3±1.3 


-19.4±1.3 


9.48±1.51 




HD 148982 


8.669±0.026 


8.381±0.036 


8.308±0.023 


-15.8±2.0 


-25.4±1.7 






HD 154922 


7.742±0.023 


7.540±0.042 


7.456±0.021 


-7.4±1.3 


-13.1±1.2 






HIP 75685 


9.186±0.024 


8.870±0.042 


8.810±0.024 


-28.8±1.5 


-19.4±1.4 


8.92±1.72 




HIP 79354 


8.090±0.023 


7.752±0.031 


7.649±0.021 


-12.1±1.0 


-29.8±1.0 


9.18±1.99 




RHS 48 


8.465±0.029 


7.830±0.053 


7.624±0.024 


-15.1±1.6 


-10.6±1.7 






TYC 0376-0769-1 


8.196±0.024 


7.649±0.031 


7.496±0.036 


-18.8±1.8 


-4.4±1.7 






TYC 0909-0125-1 


9.433±0.026 


8.987±0.022 


8.896±0.023 


-16.2±2.7 


-17.6±2.6 
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Table 2 — Continued 



Plot 


Object Name 












Parallax(3) 


ID 










[mas yr~^] 


[mas yr ^] 


[mas] 




TYC 0937-0754^1 


9.016±0.029 


8.525±0.040 


8.404±0.023 


-13.1±2.0 


-1.4±2.0 






TYC 0976-1617-1 


10.158±0.026 


9.798±0.026 


9.664±0.021 


-11.0±3.0 


-15.1±3.1 






TYC 5003-0138-1 


9.084±0.023 


8.492±0.027 


8.317±0.026 


-27.3±1.6 


-13.1±1.5 






TYC 5022-0263-1 


8.855±0.029 


8.401±0.047 


8.244±0.031 


-6.7±3.2 


-24.5±3.4 






TYC 5668-0365-1 


9.564±0.024 


9.033±0.026 


8.840±0.025 


-2.8±2.2 


-18.7±2.3 






TYC 6141-0525-1 


9.335±0.027 


8.957±0.024 


8.872±0.024 


-15.2±2.3 


-35.9±2.4 






TYC 6215-0184-1 


8.677±0.026 


8.003±0.036 


7.756±0.024 


-3.6±2.9 


-21.8±3.1 





Note. — For star #19, we adopt proper motions of its companion, #1. 

i-From 2MASS Catalog. 

^■Tycho-2 proper motions. 

^ Hipparcos parallaxes. 

* Proper motions from UCAC2. 



Table 3. Effective Temperatures and Litfiium Equivalent Widths 



Plot 


Object Name 


Li I EW 


Li I EW 


Ctmn. 


Ha 


Spectral 


Type 


TefF 


A6200/A6210 


Spectral 


Tefl 


Adopted Tgfj 


ID 




[mA] Integ. 


[mA] GFit 


[mA] 


Flag(i) 


Type 


Source(^^ 


[K] 


Line Ratio 


Type 


[K] 


[K] 


1 


TYC 6242-0104-1 


491 


487 


25 


e* 


K5 Ye 


2 


4350 


1.07 


K5 


4350 


4350 


2 


TYC 6191-0552 


481 


492 


15 


e* 


K2 


1 


4900 


1.91 


K2 


4900 


4900 


3 


TYC 6234-1287-1 


464 


452 


18 


e* 


K4 Ve 


2 


4590 


2.65 


K1.5 


4990 


4590 


4 


TYC 7312-0236-1 


434 


433 


15 


e* 


K2 Ve 


2 


4900 


2.30 


K2 


4900 


4900 


5 


TYC 7327-0689-1 


416 


414 


15 


e* 


K2 Ve 


2 


4900 


2.52 


K2 


4900 


4900 


6 


TYC 6781-0415-1 


409 


426 


13 


e* 


G9 IVe 


2 


5410 


4.69 


G9.5 


5330 


5410 


7 


TYC 6803-0897-1 


408 


413 


14 


a 








3.62 


K0.5 


5165 


5165 


8 


TYC 6214-2384-1 


397 


398 


14 


a 


Kl IV 


2 


5080 


2.48 


K2 


4900 


5080 


9 


TYC 6806-0888-1 


320 


345 


12 


a 


G8 IV 


2 


5520 


11.3 


G3 


5830 


5520 


10 


BD +04 3405B 


233 


242 


5 


a 








t 






6600* 


11 


HD 144713 


164 


193 


5 


a 


F4 


5 


6590 


t 






6590 


12 


HD 153439 


180 


198 


5 


a 


F5 V 


3 


6440 


t 






6440 


13 


HD 148396 


197 


216 


14 


a 


Kl/2 -1- F 


3 


5080 


8.09 


G8.5 


5465 


5080 


14 


TYC 7334-0429-1 


368 


378 


15 


a 


K2e 


2 


4900 


3.99 


KO 


5p50 


4900 


15 


TYC 6817-1757-1 


274 


244 


13 


e* 


KO Ve 


2 


5250 


6.03 


G9 


MIO 


5250 


16 


HD 157310 


52 


76 


1 


a 


A7 ii/ni 


5 


7850 


t 




o 

1 " 


7850 


17 


HD 142016 


20 


57 





a 


A4 IV/V 


3 


8460 


t 






8460 


18 


CD -25 11942 


307 


324 


13 


a 


KO IV 


2 


5250 


8.72 


G8.5 




5250 


19 


2MASS J17215666-2010498 


223 


226 


20 


e* 








0.58 


K7.5 


4060 


4060 


20 


TYC 6790-1227-1 


324 


344 


13 


a 


G9 IV 


2 


5410 


4.53 


G9.5 


5330 


5410 


21 


TYC 7346-1182-1 


256 


262 


12 


a 


G8 V 


2 


5520 


5.94 


G9 


5410 


5520 




BD -06 4414 






10 


c 


G5 


1 


5770 


t 






5770 




HD 137396 


47 


77 


4 


a 


F2/3 IV/V 


5 


6890 / 6740 


t 






6815 




HD 138969 


113 


109 


9 


c 


Gl V 


4 


5945 


4.68 


K9.5 


3955 


5945 




HD 140574 






1 


a 


A9 V / A3 


5/1 


7390 / 8720 


t 






8055 




HD 141612 


128 


120 


10 


a 


G5 V 


5 


5770 


4.25 


KO 


5250 


5770 




HD 141693 








a 


AO 


1 


9520 


t 






9520 




HD 141813 


201 


187 


14 


a 


G8/K2 HI + F/G 


3 


5520 / 4900 


31.4 


>G2 


>5860 


5210 




HD 142987 


150 


216 


10 


e* 


G3/6 / G5 


4/6 


5830-5700/5770 


t 






5770 




HD 143332 






5 


a 


F5 V 


5 


6440 


4.01 


KO 


5250 


6440 




HD 143358 


196 


215 


9 


a 


Gl / G2V 


3 


5945 / 5860 


t 






5900 




HD 143810 


66 


67 


6 


a 


F5/6 V 


5 


6440 / 6360 


9.25 


G8 


5520 


6400 




HD 144393 






7 


a 


F7/8 V 


4 


6280 / 6200 


13.8 


>G2 


>5860 


6240 




HD 144726 


29 


36 


5 


a 


F5 V 


5 


6440 


t 






6440 



Table 3 — Continued 



Plot Object Name 


Li I EW 


Li I EW 


Ctmn. 


Ha 


Spectral 


Type 


Tcff 


A6200/A6210 


Spectral 


Tcff 


Adopted Tcff 


ID 


[mA] Intcg. 


[mA] GFit 


[mA] 


Flag(i) 


Type 


Source(^) 


[K] 


Line Ratio 


Type 


[K] 


[K] 


HD 144732 


178 


195 


9 


a 


GO V / GO 


3/6 


6030 


27.1 


>G2 


>5860 


6030 


HD 145169 


51 


46 


10 


a 


G3 V 


5 


5830 


13.5 


>G2 


>5945 


5830 


HD 145551 


26 


25 


6 


a 


F5/6 V 


5 


6440 / 6360 


9.17 


G8.5 


5465 


6400 


HD 148982 


181 


216 


8 


a 


F8 / GO 


3 


6200 / 6030 


t 






6115 


HD 154922 




12t 


1 


a 


A7 HI 


4 


7850 


t 






7850 


HIP 75685 


99 


99 


11 


a 


G6 


1 


5700 


9.05 


G8.5 


5465 


5700 


HIP 79354 


38 


41 


8 


a 


F8 


1 


6200 


2.87 


Kl 


5080 


6200 


RHS 48 


157 


153 


15 


e* 


K2 


1 


4900 


1.77 


K2 


4900 


4900 


TYC 0376-0769-1 


12 


21 




a 








t 








TYC 0909-0125-1 


53 


54 


14 


a 








4.19 


K0.5 


5165 


5165 


TYC 0937-0754-1 


99 


97 


16 


c 








1.89 


K2.5 


4815 


4815 


TYC 0976-1617-1 






3 


a 


A9 V 


5 


7390 


t 






7390 to 


TYC 5003-0138-1 


79 


79 


16 


e*,p 








1.85 


K2.5 


4815 


4815 


TYC 5022-0263-1 








o 








t 








TYC 5668-0365-1 


61 


70 




e 








t 








TYC 6141-0525-1 








a 








t 








TYC 6215-0184-1 


49 


52 


15 


c* 


K2IVe 


2 


4900 


t 






4900 



Note. — We report here two Li A6707 measurements- "Integ." in column 3 refers to direct integration over the line profile, and "Gfit" in column 4 indicates the result 
of fitting a Gaussian to the absorption feature. In column 5, we also report contamination (denoted "Ctmn.") of the Li I line; see ^ 13.21 for description of its derivation, f 
Blended line; result indicates Gaussian feature fit to Li I in deblending. | denotes cases in which the line ratio could not be measured from the spectrum either due to extreme 
rotational broadening or the lack of prese nce of either or both l ines i n question. * Effective temperature determined via interpolation of dereddened H — K color over the 
color-effective temperature relationship of lKenvon fc HartmannI (^99^, see Fig. [2] 

^ Indicator flags are defined as follows: Absorption, a; core filling observed, c; double peaked emission, e*; P-Cygni like feature, p; emission with overlaid absorption, o. 

•^'S pectral types drawn from the following sources: typed by A.J . Weinberger using KAST low-resolution spectrograph, 1: [Torres et al.l fSACY. l2006h . 2; Michigan spectral 
atlas llHouklll98j:lHouk fc Smith-Moor3ll988l: iHouk fc Swiftlll999ll . 3, 4, and 5, respectively; HD Catalog spectral type, 6. 



Table 4. Kinematic Analysis of Membership Probability 



Plot Object Name Vsini Measured Vr Predicted Vr Ht Comovement Distance Predicted Vr /^t Comovement Distance 

ID [km s~^] [km s^^] [km s^^] [mas yr^^] Probability [pc] [km s"-*^] [mas yr^-'^] Probability [pc] 





Upper Sco 


Velocity Model 






UCL Velocity Model 






A 


HD 141569 




-6±5 


-9.4 


-4.9 ± 1.1 


15.5 


127 ± 11 


-8.1 


-0.7 ± 1.1 


95.4 


148 ± 11 


1 


TYC 6242-0104-1 


14 


-7.5±1.9 


-7.5 


-8.4± 3.6 


9.0 


228 ± 59 


-7.5 


-5.1± 3.7 


41.7 


238 ± 56 


2 


TYC 6191-0552-1 


<10 


3.3±3.0 


-6.1 


-4.2± 3.4 


56.4 


151 ± 25 


-4.7 


0.0 ± 3.4 


100.0 


171 ± 26 


3 


TYC 6234-1287-1 


22 


-5.8±1.9 


-8.7 


-1.7± 2.6 


92.0 


87 ± 9 


-8.8 


6.5 ± 2.6 


22.7 


101 ± 9 


4 


TYC 7312-0236-1 


24 


4.7±1.7 


-0.7 


-6.3± 3.3 


26.2 


140 ± 19 


1.4 


-1.8 ± 3.3 


88.7 


152 ± 20 


5 


TYC 7327-0689-1 


16 


-0.7±1.7 


-1.8 


-4.0± 3.2 


61.3 


118 ± 14 


0.0 


1.6 ± 3.2 


91.8 


130 ± 15 


6 


TYC 6781-0415-1 


29 


-3.0±2.1 


-3.4 


-2.6± 2.6 


77.8 


109 ± 11 


-1.7 


3.1 ± 2.6 


65.8 


123 ± 12 


7 


TYC 6803-0897-1 


22 


-4.2±1.9 


-4.3 


-4.3± 2.5 


46.8 


120 ± 13 


-3.3 


1.5 ± 2.5 


89.3 


134 ± 13 


8 


TYC 6214-2384-1 


15 


-3.3±1.4 


-5.5 


-6.8± 3.5 


27.8 


121 ± 17 


-4.4 


-1.2 ± 3.6 


95.6 


134 ± 18 


9 


TYC 6806-0888-1 


43 


-2.2±1.8 


-3.7 


-2.2± 3.0 


84.6 


127 ± 15 


-2.7 


3.3 ± 2.9 


64.6 


144 ± 16 


10 


BD -1-04 34056* 


39 


-14.1±2.2 


-13.9 


-0.6± 1.5 


96.2 


169 ± 23 


-14.3 


3.3 ± 1.6 


29.2^ 


207 ± 27 


11 


HD 144713 


>70 


0.6±4.4 


-4.4 


-1.9± 1.7 


72.7 


163 ± 16 


-3.1 


2.1 ± 1.7 


64.0 


184 ± 17 


12 


HD 153439 


56 


-4.3±1.6 


-5.0 


1.4± 1.6 


88.3 


130 ± 12 


-4.5 


6.8 ± 1.6 


2.0 


151 ± 13 


13 


HD 148396 


36 


-0.2±2.1 


-3.4 


-0.9± 2.4 


94.4 


231 ± 37 


-2.4 


2.1 ± 2.4 


72.6 


261 ± 41 


14 


TYC 7334-0429-1 


27 


-4.3±1.5 


-2.2 


-5.2± 2.2 


25.3 


129 ± 13 


-0.8 


0.1 ± 2.2 


99.9 


142 ± 14 


15 


TYC 6817-1757-1 


11 


8.4±1.9 


-4.8 


-7.4± 2.8 


3.8 


379 ± 99 


-4.1 


-5.4 ± 2.7 


16.5 


382 ± 90 


16 


HD 157310 


>70 


-22.4±2.3t 


-13.9 


-1.6± 1.5 


74.1 


204 ± 28 


-14.3 


1.7 ± 1.5 


67.7 


245 ± 32 


17 


HD 142016 




-17.8±2.3t 


-2.3 


-6.2± 1.2 


22.6 


85 ± 7 


-0.7 


1.7 ± 1.2 


87.3 


94 ± 6 


18 


CD-25 11942 


53 


-6.5±1.5 


-5.7 


-2.1± 2.0 


78.0 


127 ± 12 


-5.4 


3.6 ± 2.0 


43.4 


144 ± 13 


19 


2MASS J17215666-2010498 


<10 


-7.3±1.9 


-7.5 


-8.4± 3.6 


9.0 


228 ± 59 


-7.5 


-5.1 ± 3.7 


41.7 


238 ± 56 


20 


TYC 6790-1227-1 


33 


2.0±2.0 


-2.6 


-6.1± 2.7 


23.7 


121 ± 13 


-1.0 


-0.7 ± 2.6 


98.1 


133 ± 13 


21 


TYC 7346-1182-1 


28 


3.7±1.7 


-3.6 


-4.8± 2.3 


31.8 


129 ± 13 


-2.8 


0.7 ± 2.3 


97.1 


143 ± 14 



Note. — fVery broadened, featureless spectra. Radial velocity derived via centroid measurement of the Ha line. Uncertainty reflects an assumed centroid measurement error of 0.05A. For 
a full discussion of error analysis, see ^ 13.31 
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20 340 320 300 280 260 

in 



Fig. 1. — A galactic coordinate map of HD 141569 and nearby associations is presented. 
Open circles indicate our sample stars, while dashed boxes indicate regions studied by 
de Zeeuw et al.l (1l999l ). Det ails of objects plotted herein are to be found in the foUow- 



de Zeeuw et al. 



ing papers: Upper S corpius (iPreibisch et 3.1.1120021) . UCL and L CC (iMamaiek et al 



Zuckerman et al 



99^. TW Hya jMamaiekl boosh. n a nd e Cha JZuckerman fc Song 



2002 



2004 



200l|) and IC 2602 flRobichon et al.l 119991 ). In all subsequent plots, the 



letter A denotes the position of HD 141569. 



-24- 



8000 



7000 



6000 



0) 



5000 



4000 



3000 



T — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I 1 — I — I — I — I — I — I — I — I — I — I — T 

o Observed target ^ 
▲ Spectral Type Standard 

O 



~i — I — I — I — I — I — I — I — I — I — I — I — I — I — r 



o 



O O I 
g O ©I 



O 



CD 



O 



UU 1 

O p 



o (ffio Q OA ^' 



OA 



o / 

/ 

/ 



_] I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I L 



0.4 



0.3 



0.2 



0.1 



0.0 



-0.1 



(H-K) 



2MASS 



Fig. 2. — Color- Tgfj relationship from lKenyon fc Hartmanru (119951 ) for 2 MASS (H—K) colors. 
The s olid line represents the main-sequence relationship as defined in iKenyon fc Hartmann 
(119951 ). Open circles represent literature spectral types for our sample as reported in column 6 
of Table [31 Filled triangles represent literature spectral types for the standard stars observed 
by us (Table [1]). 
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Fig. 3. — Near-infrared JHK colors from 2MASS. For visual clarity, we display here only 
the 21 stars that we identify as bein g lithium rich (see ^3.21 and Fig. HI). Solid lines represent 
the dwarf and giant sequences from Bessell ( 1991 ) and lSessell fc Brett ( 1988 ). Dereddened 
objects are plotted with observed and dereddened color s connected by a thin , solid line 
parallel to a reddening vector defined by = 1.95 (IBessell fc BrettI Il988l ). Numbers 

next to the points are provided for ease in identifying the objects in the data tables. Stars 
with Ha in emission are indicated with an X over the point. 
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Fig. 4. — Fe contaminat ion corrected L i I EW as a function of Tog, plotted with the IC 
2602/IC 2391 data from iRandich et al.l (119971 . 12OOII ). Stars possessing enough hthium to 
be above the selection threshold are likely to be pre-main-sequence, and are designated 
the "High Li" sample. Stars below the threshold but above the horizontal line at 200mA 
constitute the "Moderate Li" sample; these stars have Li EWs in the upper envelope of 
the IC2602/IC2391 locus. We include objects 16 and 17 in the moderate sample to be 
conservative as Li I is not a good age indicator in higher mass stars. Numbers next to the 
points are provided for ease in identifying the objects in the data tables. In two cases we 
show both the measured EW as well as the doubled value (indicated by a "b" after the 
number identifier), these objects may be suffering Li line filling (see §3.2). Stars showing Ha 
in emission are indicated with an X over the point. 
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Fig. 5. — A cropped area about HD 141569 plotted with high and moderate hthium sample. 
Symbols are defined as in Fig. [H 
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Predicted Upper Sco Radial Velocity 




30 20 10 350 340 330 

Fig. 6. — Radial velocities for the Lithium selected sample plotted as a function of galactic 
longitude. Also shown for comparison, the radial velocities predicted from the Upper Sco 
velocity vector at a galactic latitude of +20°are projected in galactic longitude. Shown are 
la uncertainties; all high Lithium stars are consistent in radial velocity space with Upper 
Sco within 3a. 




Fig. 7. — HR Diagram for sample stars with the pre-main-sequence tracks of iBaraffe et al. 
(Il998l ) overplotted. Isochrones shown are (fro m top to bottom) 1, 3, 10, 20 , 30, and 100 
Myr. Also shown is the 100 Myr isochrone of iD'Antona fc Mazzitellil (119971 ) extending to 
higher masses. This isochrone and the 10 Myr isochrone are highlighted for reference. The 
upper panel shows the stars using M^s derived using the d istance to HD 141569 from the 



second Hipparcos data release, 116±8 pc (Ivan Leeuwen 



200 



3) ■ In the mid dle pane, we apply 



the mean distance to Upper Scorpius, 145±2 pc (jde Zeeuw et al.l 119991 ). The lower panel 
was generated using distances derived from the Upper Sco velocity model (see text). The 
sample's isochronal age appears consistent with the high lithium abundances which indicate 
ages <30 Myr (see Fig. H]). 
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Table 5. Association Notes 



Plot 


Object 


Spatial 


Velocity Model / 


US/UCL Vr 


Determined 


ID 


Name 


Match 


Probability 


Prediction 


Membership 


1 


TYC 6242-0104-1 


? 


US / 9% 


la/la 


I 


2 


TYC 6191-0552 


US 


US / 56% 


no/3cr 


US 


3 


TYC 6234-1287-1 


? 


US / 92% 


2a /2a 


US? 


4 


TYC 7312-0236-1 


UCL 


UCL / 89% 


no/2(T 


UCL 


5 


TYC 7327-0689-1 


UCL 


UCL / 92% 


la/la 


UCL 


6 


TYC 6781-0415-1 


US 


US / 78% 


la/la 


US 


7 


TYC 6803-0897-1 


us 


US / 47% 


la/la 


US 


8 


TYC 6214-2384-1 


us 


US / 28% 


2a/la 


US 


9 


TYC 6806-0888-1 


us 


US / 85% 


la/la 


US 


10 


BD +04 3405B 


7 


US / 96% 


la/la 


I 


11 


HD 144713 


us 


US / 73% 


2a/la 


us 


12 


HD 153439 


7 


US / 88% 


la/la 


us? 


13 


HD 148396 


us 


US / 94% 


2a/2a 


us 


14 


TYC 7334-0429-1 


UCL 


UCL / 100% 


2a /3a 


UCL 


15 


TYC 6817-1757-1 


us 


US / 4% 


no/ no 


I 


16 


HD 157310 


7 


US / 78% 


no/ no 


I 


17 


HD 142016 


UCL 


UCL / 87% 


no/ no 


UCL 


18 


CD-25 11942 


7 


US / 78% 


2a/2a 


us? 


19 


2MASS J17215666-2010498 


7 


US / 9% 


la/la 


I 


20 


TYC 6790-1227-1 


UCL 


UCL / 98% 


Za/2a 


UCL 


21 


TYC 7346-1182-1 


US 


US / 32% 


no/ no 


US 



Note. — In column three, we note with which moving group each object is spatially consistent. Column 
four summarizes the results of velocity vector modeling and shows the vector with which each object 
had the highest comovement probability. In column five, we report the models which predicted radial 
velocities within 2a of our measured values, and in the rightmost column we comment on membership; 
I denotes "indeterminate". See text (§[5]l for further discussion of our final membership determinations. 
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